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38 Abstract In this study, we assess changes of aerosol optical depth (AOD) and direct 

39 radiative forcing (DRF) in response to the reduction of anthropogenic emissions in four 

40 major pollution regions in the northern hemisphere by using results from 10 global 

41 chemical transport models in the framework of the Hemispheric Transport of Air 

42 Pollution (HTAP). The multi-model results show that on average, a 20% reduction of 

43 anthropogenic emissions in North America, Europe, East Asia and South Asia lowers the 

44 global mean AOD and DRF by about 9%, 4%, and 10% for sulfate, organic matter, and 

45 black carbon aerosol, respectively. The impacts of the regional emission reductions on 

46 AOD and DRF extend well beyond the source regions because of intercontinental 

47 transport. On an annual basis, intercontinental transport accounts for 10-30% of the 

48 overall AOD and DRF in a receptor region, with domestic emissions accounting for the 

49 remainder, depending on regions and species. While South Asia is most influenced by 

50 import of sulfate aerosol from Europe, North America is most influenced by import of 

5 1 black carbon from East Asia. Results show a large spread among models, highlighting the 

52 need to improve aerosol processes in models and evaluate and constrain models with 

53 observations. 

54 


2 



55 1 . Introduction 

56 Anthropogenic aerosols make significant contributions to the global mean radiative 

57 forcing (RF) of climate [Forster et al., 2007] by scattering and absorbing solar radiation 

58 (so-called “aerosol direct effects”) [McCormick and Ludwig, 1967] and modifying cloud 

59 properties, amount, and evolution (collectively referred to as “aerosol indirect effects”) 

60 [Twomey, 1977; Gunn and Phillips, 1957; Albrecht, 1989] . RF is a measure of the 

61 change of net radiation (incoming minus outgoing) at the top of atmosphere (TOA), at the 

62 surface, or within the atmosphere, due to perturbations in atmospheric compositions or 

63 surface properties. On a global average basis, the sum of direct and indirect RF at TOA 

64 by anthropogenic aerosol is estimated to be -1.2 W m" 2 [-2.4 to -0.6 W m" 2 ] (cooling) 

65 over the period of 1750-2000, which is significant compared to the positive (warming) 

66 forcing of +2.63 [+0.26] W m" 2 by anthropogenic long-lived greenhouse gases over the 

67 same period [Forster et al., 2007], In heavily polluted regions, aerosol cooling 

68 overwhelms greenhouse warming [Ramanathan et al, 2001; Li et al., 2010] . At the 

69 surface, aerosol RF can be much stronger than that at TOA because of aerosol absorption 

70 [Satheesh and Ramanathan, 2000] . Currently, uncertainties associated with aerosol RF 

7 1 make the largest contribution to the overall uncertainty in anthropogenic radiative forcing 

72 of climate [Forster et al., 2007] . Because of the significant role of aerosols in modulating 

73 the Earth’s radiative budget, it is necessary from both scientific and policy perspective to 

74 assess how emission changes associated with economic development and 

75 regional/national regulations will influence the aerosol radiative forcing. 

76 
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77 The response of global aerosol RF to a fractional change of anthropogenic emissions 

78 would depend on the source locations, magnitude, and composition of emitted aerosols 

79 and aerosol precursors. While scattering aerosols like sulfate cause a cooling effect, 

80 strongly absorbing black carbon aerosols cause warming. Aerosol RF is also determined 

81 by several environmental factors such as surface albedo and meteorological conditions 

82 (in particular the amount and distribution of clouds and winds). With these differences, 

83 the aerosol RF per unit change of regional emissions would differ from one region to 

84 another. It is also important to note that the impact of a regional emission reduction is not 

85 confined to the region itself. Instead regional emission reductions can have far reaching 

86 impacts on RF in downwind regions, because of intercontinental transport of 

87 anthropogenic aerosols. Long-range transport has been observed by long-term surface 

88 monitoring networks [Prospero et al., 2003; VanCuren, 2003; Fischer et al., 2010], in- 

89 situ measurements from intensive field campaigns [Ramanathan et al., 2007; Clarke and 

90 Kapustin, 2010] , and satellite observations [Yu et al., 2008; Rudich et al., 2008; Dirksen 

91 et al., 2009] backed by model simulations [Heald et al, 2006; Chin et al., 2007; Hadley 

92 et al, 2007] . Thus how a region is influenced by extraregional emissions is of particular 

93 importance for formulating an effective strategy for mitigating regional climate change 

94 and combating air pollution. The United Nations Economic Commission for Europe has 

95 established a Task Force on Hemispheric Transport of Air Pollutants (HTAP) to 

96 understand the growing body of scientific evidence of intercontinental transport and 

97 assess its impacts on air quality, climate, and ecosystems [HTAP, 2010] . 

98 
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99 Modeling studies can offer valuable insights into the relative significance of aerosols 

100 from different regions in influencing climate and important implications for formulating 

101 effective emission control strategies. Several recent studies have assessed how aerosols 

102 emitted in a region or from specific sectors could affect climate in downwind regions 

103 [e.g., Reddy and Boucher, 2007; Koch et al., 2007; Shindell et ah, 2008a, 2008b], Other 

104 studies have shown that large intennodel differences exist in aerosol lifetime cycle and 

105 radiative effect [Kinne et al., 2006; Schulz et al., 2006; Textor et al., 2006], which might 

106 undennine the robustness of the results from single model or very limited number of 

107 models. Complementary to the prior studies, we use in this study an ensemble of 10 

108 global chemical transport or general circulation models to assess the change of global and 

109 regional aerosol optical depth (AOD) and direct radiative forcing (DRF) in response to 

110 20% reductions of emissions from four major polluted regions in the Northern 

1 1 1 Hemisphere. These models conducted aerosol source-receptor (S/R) relationship 

1 12 experiments in the HTAP coordinated studies (HTAP 2010). These multi-model S/R 

113 experiments allow us to examine a probable range of contributions of intercontinental 

114 transport (ICT) relative to intra-regional emissions in determining regional AOD and 

115 DRF and help characterize the robustness of the results. Fry et al. [2012] conduct similar 

116 analysis of ozone radiative forcing due to 20% reductions of ozone precursor emissions 

117 using results from multiple HTAP models. 

118 

119 The rest of this paper is organized as follows. Section 2 describes the S/R model 

120 simulations and analysis methodology, including AOD from the HTAP anthropogenic 

121 S/R experiment used in this study and an estimate of the aerosol DRF. Section 3 presents 
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122 a multi -model analysis of the impacts of 20% reductions of regional anthropogenic 

123 emissions on global and regional AOD and DRF and the role of ICT relative to intra- 

124 regional emissions. This assessment does not address aerosol indirect effects, or warming 

125 effects resulting from BC deposition on snow and ice. We also neglect interactions of 

126 aerosols with thermal infrared radiation, as anthropogenic aerosols have relatively small 

127 size and their interactions with infrared radiation are minor. Finally, we neglect the 

128 radiative forcing of gas-phase components that may have been influenced by these 

129 emissions, which were modeled by Fry et al. (2012). Major conclusions from the analysis 

130 are summarized and discussed in Section 4. 

131 

132 2. Description of model simulations and analysis methodology 

133 2.1. HTAP modeling experiments on S/R relationships 

134 We use the results from an ensemble of 10 global chemical transport or general 

135 circulation models, listed in Table 1, to evaluate changes in AOD and DRF in response 

136 to a 20% reduction of anthropogenic emissions in four major pollution regions in the 

137 Northern Hemisphere under the framework of the HTAP S/R modeling experiments 

138 [Fiore et al., 2009; HTAP, 2010] . The S/R experiments include a baseline simulation and 

139 four perturbation simulations, for which each model submitted monthly AOD fields. For 

140 each model, the baseline simulation (SRI) is conducted using emissions and meteorology 

141 for 200 1 . Note that individual models use their own preferred anthropogenic and natural 

142 emissions. Each of the four perturbation runs (SR6) represents a 20% reduction in 

143 anthropogenic emissions (including biomass burning and n=both gas-phase and aerosol 

144 components) in one of four major pollution regions in Northern Hemisphere, namely 
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145 North America (NA), Europe (EU), East Asia (EA), and South Asia (SA), as illustrated in 

146 Figure 1. These perturbation model experiments are denoted as SR6NA, SR6EU, 

147 SR6EA, and SR6SA, respectively. Also shown in Figure 1 are regional and annual 

148 anthropogenic emissions of sulfur dioxide (SO 2 ), primary particulate organic matter 

149 (POM), and black carbon (BC) from 8 models (emissions for ECHAM5 and HadGEM2 

150 were not archived and can not be retrieved for this analysis). Clearly anthropogenic 

151 emissions show large regional differences. For example, South Asia has the least SO? 

152 emissions that are a factor of 2-4 smaller than the other regions. East Asia has the largest 

153 BC emissions that are nearly double the emissions in Europe or South Asia and more 

154 than triple the emissions in North America. The regional differences in emissions in 

155 combination with the proximity of regions (e.g., the SA region is adjacent to and 

156 downwind of the EU region), ICT pathways, and removal/production mechanisms during 

157 the transport, would detennine the relative roles of ICT and regional emissions, as we 

158 will discuss later in the paper. 

159 

160 2.2 AOD from HTAP S/R experiments 

161 The major HTAP model outputs used in this study are monthly average AOD (x) at 550 

162 mn for sulfate (SO4), POM, and BC from individual models. Figure 2 shows annual 

163 average AOD in ambient conditions for a combination of SO4, POM, and BC from 9 

164 models. ECHAM5 simulations are not shown here because the model calculates AOD for 

165 dry sulfate, POM, and BC. Although the model also provides water optical depth 

166 associated with total aerosol (e.g., due to aerosol humidification), it is impossible, 

167 without further uncertain assumptions, to partition the water optical depth into that 
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168 associated with individual aerosol components. Note that these AOD outputs include both 

169 anthropogenic and natural contributions (e.g., DMS derived sulfate AOD and biomass 

170 burning AOD, among others). The AOD distributions clearly show several hotspots 

171 representing well-known industrial pollution regions (e.g., East and South Asia, Western 

172 Europe, and eastern US) and biomass burning regions (e.g., equatorial Africa in DJF, 

173 South America and southern Africa in JJA and SON). In Table 2, the global annual mean 

174 AOD is 0.0335 ± 0.0134 (average ± standard deviation of the 9 models), 0.0108 ± 

175 0.0047, and 0.0022 ± 0.0024 for SO4, POM, and BC, respectively. Clearly AOD shows 

176 large model diversity with the standard deviation equivalent to about half of the multi - 

177 model average. These global annual means are lower than those from multiple model 

178 simulations with harmonized emissions under the framework of the Aerosol Comparisons 

179 between Observations and Models (AeroCom) [Schulz et a/., 2006] by 10%, 28%, and 

180 8% for sulfate, POM, and BC, respectively. 

181 

182 Figure 3(a) shows the 9-model median AOD for sulfate, POM, and BC combined as 

183 derived from the SRI baseline simulations and its comparisons with anthropogenic AOD 

184 derived from MODIS over-ocean measurements using two different methods described 

185 by Yu et al. [2009] (Figure 3b, denoted as MODIS-YU09) and Bellouin et al. [2008] 

186 (Figure 3c, denoted as MODIS-BE08). MODIS-YU09 anthropogenic AOD is derived 

187 from MODIS over-ocean retrievals of AOD and fine-mode fraction (FMF) using the 

188 representative FMF values for individual aerosol types, which are determined from 

189 MODIS observations in selected regions where the specific aerosol type predominates 

190 [Yu et al., 2009] . MODIS-BE08 is derived from the same MODIS measurements, but 
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191 with the use of prescribed, in situ measurement-based thresholds of FMF for different 

192 aerosol types (which are different from those derived from MODIS observations as in Yu 

193 et ah, 2009) aided by the satellite observed absorbing aerosol index to separate 

194 anthropogenic aerosol from mineral dust and sea salt [Bellouin et al., 2008]. Deriving 

195 anthropogenic AOD from satellite measurements over land is not currently feasible. As 

196 shown in Figure 3, anthropogenic AOD of MODIS-YU09 is generally larger than 

197 MODIS-BE08 AOD in northern hemispheric mid-latitudes. In the tropics, MODIS-BE08 

198 AOD is somewhat higher than MODIS-YU09. Except for trans-Pacific transport in 

199 spring, the cross-ocean transport from major industrial pollution and biomass burning 

200 regions is generally more extensive in the HTAP model simulations than both MODIS- 

201 based estimates. In particular, the models simulated significant cross- Atlantic transport of 

202 tropical African smoke to South America in boreal winter, which is not clearly seen in the 

203 MODIS-based estimates. Figure 4 further compares zonal variations of HTAP seasonal 

204 AOD for the sulfate, POM and BC mixture against that of the MODIS anthropogenic 

205 AOD over the ocean in mid-latitudes (20-60°N) of the Northern Hemisphere, where the 

206 major intercontinental transport paths occur. Significant regional and seasonal differences 

207 exist between the HTAP models and MODIS observations and between MODIS-YU09 

208 and MODIS-BE08. In the North Atlantic, the MODIS anthropogenic AOD from both 

209 methods is often near the low bound of 9 HTAP models. In the North Pacific, MODIS- 

210 YU09 anthropogenic AOD generally agrees well with the median of HTAP models. On 

211 the other hand, MODIS-BE08 is consistently lower than MODIS-YU09 and generally 

212 entails the low bound of the HTAP model simulations, particularly in the northeastern 

213 Pacific. The differences between the two MODIS-based estimates of anthropogenic AOD 
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214 result from the differences in specifics of the approaches, as discussed in Yu et al. 

215 [2009] . The pronounced differences between the two MODIS-based estimates of 

216 anthropogenic AOD manifest the urgent need to develop independent approaches based 

217 on other satellite measurements, and more importantly, to conduct in situ measurements 

218 of anthropogenic AOD over large scales [Yu et al., 2009] . 

219 

220 2.3 Estimate of the aerosol direct radiative forcing 

221 Because the aerosol direct radiative forcing or the aerosol optical properties (i.e., single- 

222 scattering albedo and asymmetry factor) for calculating the forcing are not archived under 

223 the HTAP S/R experiment protocol, we estimate here the aerosol DRF for each model 

224 and component (i.e., sulfate, POM, or BC) as follows: 

225 DRFi(x,y,t) = AOD l( x,y,t) •NDRF i (x,y,t) (i = S0 4 , POM, and BC) (1) 

226 where x, y, and t represents longitude, latitude, and month, respectively, and NDRF is 

227 the normalized DRF with respect to AOD at 550 nm [Zhou et al., 2005] . In this study we 

228 derive monthly average NDRF; over each model grid cell by dividing the DRF by the 

229 AOD, using 2001 monthly average AOD and direct radiative forcing calculated from the 

230 Goddard Chemistry Aerosol Radiation and Transport (GOCART) model [Chin et al., 

23 1 2002] . We then apply this NDRF; to the monthly average component AOD; from the 

232 other models in the HTAP S/R experiments. We also estimate DRF for an external 

233 mixture of SO4, POM and BC by simply summing up the DRF for individual 

234 components. Note that RF for an internal mixture could differ significantly from that for 

235 an external mixture [Jacobson, 2001; Chung and Seinfeld, 2002] . Although DRF does 

236 not increase with AOD in a fully linear manner over the whole range of AOD, the use of 
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237 the above linear formula to derive DRF and its change in response to a 20% reduction of 

238 emissions would not introduce large uncertainties with respect to regional differences in 

239 DRF [Zhou et a/., 2005; Anderson et al, 2005] . 

240 

24 1 The GOCART model currently prescribes particle size distributions and refractive indices 

242 for individual components based on the Optical Properties of Aerosols and Clouds 

243 (OP AC) database [Hess et al., 1998] . Aerosol properties such as AOD, single scattering 

244 albedo, and asymmetry factor are then calculated using the Mie code [Chin et al., 2002, 

245 2009] . These aerosol optical properties along with surface albedos and cloud fields from 

246 the Goddard Earth Observing System (GEOS) Data Assimilation System - Version 4 are 

247 then used to drive the NASA Goddard radiative transfer model [Chou et al., 1998], The 

248 time step for the radiative transfer calculations is 30 min, which can adequately capture 

249 the dependence of DRF on solar zenith angle [Yu et al., 2004] . For a specific component 

250 (e.g., SO 4 , POM, or BC), DRF is calculated as the difference of net downward radiative 

25 1 flux between a radiative transfer calculation including all aerosol components and one 

252 with the specific component excluded. This DRF is therefore different from the DRF 

253 estimated in Forster et al. [2007] , where the reference is pre-industrial aerosols. In this 

254 study, the GOCART DRF is calculated for solar radiation only and averaged over a 24- 

255 hour period. GOCART calculations of DRF have been evaluated against remote sensing 

256 measurements and other model simulations [e.g., Yu et al., 2004; Yu et al., 2006], 

257 

258 Table 3 lists the annual mean NDRF in the four defined regions and globally. For purely 

259 scattering SO4, NDRF is nearly the same at the top of atmosphere (TOA) and at the 
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260 surface. Because POM is partially absorbing in the UV range [Chin et ah, 2009], the 

261 surface NDRF is more negative than the TOA NDRF. BC aerosol is strongly absorbing 

262 over the whole solar spectrum, thus its TOA forcing is positive (i.e., TOA wanning), 

263 while the surface forcing is strongly negative (i.e., surface cooling). Clearly BC aerosol is 

264 much more effective in interacting with solar radiation than SO4 and POM, although it is 

265 typically associated with lower AODs. For a specific component, the difference in NDRF 

266 among regions is generally within 20-30%, due to the combined effects of differences in 

267 solar zenith angle, surface albedo, and cloud fields [Yu et al., 2006] . The global annual 

268 mean all-sky TOA NDRF is -24.2, -30.0, and +85.9 WmV for S0 4 , POM, and BC, 

269 respectively. These values fall within the ranges reported in the literature, i.e., -10 ~ -32 

270 Wm'V 1 for S0 4 , -5 ~ -38 WmV for POM, and +22 ~ +216 Win'Y 1 for BC [CCSP, 

271 2009]. 

272 

273 To apply the GOCART-based monthly NDRF; to the HTAP models, we regrid the 

274 component AOD; from each model to a unified spatial resolution of 2°x2.5° and use eq. 

275 (1) to estimate the corresponding DRF,. For the SRI experiment and on a global and 

276 annual basis, the 9 models give the mean all-sky TOA DRF of -0.67 ± 0.28, -0.26 ± 0.13, 

277 and +0.22 ±0.13 Wm" 2 for S0 4 , POM, and BC, respectively. Correspondingly, the 

278 respective DRF at the surface and in all-sky (including both cloudy and cloud-free) 

279 condition is -0.64 ± 0.26, -0.37 ± 0.19, and -0.43 ± 0.22 Wm" 2 . When the anthropogenic 

280 fraction of AOD, e.g., 55%, 53%, and 100% respectively for S0 4 , POM, and BC on a 

28 1 global and annual average basis [Schulz et al., 2006] is considered, our DRF estimates 

282 shown above are in good agreement with anthropogenic TOA DRF of -0.35 ±0.15 Wm" 2 
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283 (SO 4 ), -0. 13 ± 0.05 Wm" 2 (POM), and +0.25 ± 0.09 (BC) Wm' 2 reported by Schulz et al. 

284 [2006] . The aerosol direct radiative forcing exhibits large regional and seasonal 

285 variations, as shown in Figure 5 for the all-sky DRF and for the external mixture of SO4, 

286 POM, and BC, as derived from the 9-model SRI simulations. Over major industrial 

287 pollution regions in northern hemispheric mid-latitudes, aerosol DRF can be as large as - 

288 6 Wm" 2 at TOA and -25 Wm 2 at the surface. The forcing is greater in JJA and MAM 

289 than in DJF and SON. Similar magnitudes of DRF are seen in the tropics where biomass 

290 burning smoke is dominant. The seasonal variations of DRF reflect the peak seasons of 

291 biomass burning, e.g., DJF in the Sahel and JJA and SON in southern Africa and South 

292 America. 

293 

294 We note that the use of GOCART-based NDRF is likely to understate the model diversity 

295 in DRF, in comparison to that derived from full radiative transfer calculations by 

296 participating models. Although other aerosol models currently use the OP AC to 

297 characterize aerosol microphysical and optical properties just like GOCART, models can 

298 differ substantially in aerosol vertical distributions, meteorological fields (such as relative 

299 humidity, cloud distributions, and surface albedos), and radiative transfer schemes. These 

300 differences combined may lead to large model differences in DRF. This could 

301 particularly be the case for all-sky TOA DRF by BC, which depends strongly on the 

302 vertical distributions of aerosol and clouds, as evidenced by the reported wide range of 

303 DRF [CCSP, 2009]. 

304 
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305 3. Influences of regional anthropogenic emission reductions on AOD 

306 and DRF 

307 In this section we present how the 20% reduction of regional anthropogenic emissions 

308 changes AOD and DRF by analyzing differences between a set of SR6 experiments 

309 (SR6NA, SR6EU, SR6EA, and SR6SA) and the SRI experiment. We examine first the 

310 changes in global mean AOD and DRF and then the spatial extents of AOD and DRF 

311 changes resulting from regional anthropogenic emission reductions. The relative roles of 

312 ICT and regional emissions are assessed with a concept of the relative annual 

313 intercontinental response (RAIR) (as defined in HTAP, 2010). For the reasons discussed 

3 14 earlier in 2.2, we exclude ECHAM5 results in the analysis of AOD and DRF magnitude. 

315 Because RAIR is less sensitive to water optical depth than AOD is, we include ECHAM5 

316 results in the RAIR analysis. 

317 

318 3.1. Global mean 

319 Just as aerosol radiative forcing varies regionally, the effects of changes in emissions of 

320 aerosols and their precursors on global mean DRF also varies with region. Tables 4 and 

321 5 list respectively changes of global annual average AOD and all-sky DRF in response to 

322 the 20% reduction of anthropogenic emissions in the four HTAP anthropogenic source 

323 regions. The combined impact of the 20% reduction of emissions in all of the four 

324 regions is to decrease AOD and DRF by about 9%, 4%, and 10% for sulfate, POM, and 

325 BC, respectively. Relative contributions from individual regions vary considerably. For 

326 sulfate, the change of global average AOD and DRF due to the reduction of SO 2 

327 emissions in South Asia is substantially smaller than that due to the emission reductions 
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328 from the other regions. This is mainly because SO 2 emissions in the SA region are much 

329 smaller than those in other regions (see Figure 1). For POM, the reductions of global 

330 mean AOD and DRF due to the regional emission reduction are generally consistent with 

33 1 regional differences in primary POM emissions. For BC, the reduction of emissions in 

332 East Asia makes the largest contribution to the change of global average DRF, mainly 

333 because of the highest BC emissions among the four regions. For the external mixture of 

334 SO4, POM, and BC, the 20% reductions in the anthropogenic emissions in the four 

335 regions collectively yield the respective reductions of 0.0035, 45.5 mWm , (1 mWm' 2 = 

336 10" 3 Wm" 2 ) and 104. 1 mWm' 2 for global mean AOD, all-sky DRF at TOA and at the 

337 surface, which represents respectively about 8%, 6% and 7% reduction from the baseline 

338 simulation. Note that uncertainties are larger than given by the standard deviation in 

339 Table 5, due to the use of NDRF from a single model, as discussed earlier. 

340 

341 Will the same amount of emission reduction in different regions introduce the same 

342 change in global mean aerosol radiative forcing? To address this question, we calculate 

343 the forcing efficiency with respect to emissions by normalizing the global annual mean 

344 DRF difference between the SR6 and SRI runs by the change (20%) in regional 

345 emissions for each model, as shown in Table 6. For SO4 and POM, we use the emissions 

346 for SO 2 and primary POM respectively to calculate the forcing efficiency. These derived 

347 forcing efficiency numbers should be considered as rough estimates, because some 

348 models with fully coupled chemistry include changes of SO4 resulting from the 

349 reductions of other anthropogenic emissions (Fry et ah, 2012) and a fraction of POM is 

350 secondary aerosol formed from volatile organic carbon emissions. Clearly the forcing 
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351 efficiency for BC is an order of magnitude higher than that for SO4 and POM, suggesting 

352 that BC emission control is a candidate for producing a large, near-term cooling that 

353 counters the warming of greenhouse gases. Despite wide model differences in the 

354 absolute value of forcing efficiency, several regional dependences appear to be rather 

355 robust among the models. As shown in Figure 6, all eight models consistently yield a 

356 smaller sulfate forcing efficiency for EA emissions than for EU emissions, which is 

357 consistent with the oxidant limitation in EA that reduces the efficiency of SCE-to-sulfate 

358 transformation [Koch et al., 2007] . The models also consistently give the largest BC 

359 forcing efficiency for EU emissions. These results may have important implications for 

360 regional emission controls and their influences on global climate, which will be discussed 

361 in section 4. 

362 

363 3.2. Spatial extents and the role of ICT 

364 The spatial extents of the AOD and DRF due to the reduction of regional emissions and 

365 the role of aerosol ICT are revealed by differencing SR6 and SRI experiments. Figure 7 

366 shows the 9-model annual average changes of AOD and all-sky DRF by the external 

367 mixture of SO 4 , POM, and BC, resulting from the 20% reduction of anthropogenic 

368 emissions over the four regions. Clearly emissions from North America, Europe, and East 

369 Asia exert significant DRF on intercontinental and even hemispheric scales. Emissions 

370 from South Asia have relatively small impacts mainly over the Indian Ocean and the 

371 tropical Pacific. The influences of regional emission reduction on AOD and DRF also 

372 depend on season. As an example, Figure 8 shows seasonal variations of changes of the 

373 all-sky surface DRF between SRI and SR6. Similar seasonal variations occur for AOD, 
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all-sky DRF at TOA and clear-sky DRF, which will not be discussed separately. For 
North America, Europe, and East Asia, the 20% reduction of emissions in each region 
decreases the direct forcing (i.e., less negative) by a larger amount and over more 
extensive areas in summer and spring than in winter and fall. For comparison, the 
influence of the 20% reduction of emissions in South Asia shows relatively small 
seasonal variations. 

To quantify the role of aerosol ICT in affecting regional climate forcing, we adopt a 
concept of the RAIR as defined in HTAP [2010] . For AOD, RAIR in a receptor region i 
is expressed as follows: 


RAIRi = 


ZjSAODjj 
SAOD a+Y,j SAODji 


(j*0 


where index j represents a source region. 6AOD,, represents a change of AOD in the 
receptor/domestic region i due to the emission reduction in the region itself and hAOD,, 
represents a change of AOD in the receptor i induced by the emission reduction in a 
source region j outside of the receptor (or foreign region). Similarly, RAIR can be 
defined for DRF, near-surface concentration, and surface deposition. By definition, RAIR 
in a receptor region represents the percentage contribution of the intercontinental 
transport of foreign emissions relative to the sum of foreign and domestic emissions. 

We calculate the RAIR for AOD and DRF for individual models and then obtain multi- 
model averages, as shown respectively in Figure 9 and 10 with the columns representing 
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397 10-model (including ECHAM5) average and error bars indicating ±1 standard deviation 

398 of the model simulations. Clearly values of RAIR depend on both region and component. 

399 For all regions and components, the import from intercontinental transport is significant 

400 but local emissions remain the main contributor, with RAIR ranging from 10 to 30%. 

401 South Asia is most influenced by the import of sulfate aerosol, and North America is 

402 most influenced by the import of BC, followed by POM. For sulfate, POM, and BC 

403 combined, South Asia DRF is most strongly influenced by foreign sources, and North 

404 America and Europe is the least. These rankings also reflect the strength of local 

405 emissions relative to world emissions, as discussed in section 2. Interestingly, RAIR 

406 values for AOD and DRF of sulfate (10-30%) are consistently smaller than that for 

407 aerosol column loading (i.e., 24-37%) [see Table 4.4 in HTAP 2010], For POM, RAIR 

408 values (16-18%) for AOD and DRF in NA and EU are also somewhat smaller than the 

409 corresponding RAIR values (21-23%) for the column mass loading. These differences 

410 probably stem from transported sulfate and POM aerosols experiencing lower relative 

411 humidity at high altitude (resulting in lower optical depth) than local aerosols, which 

412 generally remain at lower altitudes. Finally, model variability is very large, highlighting 

413 the significant uncertainties in modeling aerosol processes. 

414 

415 The contributions of emission reductions in a receptor region depend on source regions. 

416 Table 7 lists the response of annual mean AOD and all-sky DRF in receptor regions to a 

417 20% reduction of anthropogenic emissions in the four source regions. Here we combine 

418 SO4 and POM together, in which SO4 makes a dominant contribution of 69-94%, 

419 depending on regions. BC is listed separately from the SO4 and POM combined, because 
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420 strongly absorbing BC has different effects on solar radiation than do weakly absorbing 

42 1 sulfate and organic matter. Not surprisingly, the influences on AOD and DRF from the 

422 reduction of domestic emissions is about an order of magnitude larger than from any 

423 foreign source region. Foreign source regions also differ in contributing to the total 

424 import for a specific receptor region. The intercontinental transport of EA emissions 

425 accounts for the largest fraction of total import in NA, i.e., 60% and 70% for SO4+POM, 

426 and BC, respectively. The NA emissions make the largest contribution (57%) to the total 

427 import of SO4+POM into EU. For BC aerosols imported to EU, on the other hand, the EA 

428 contribution of 45% exceeds the NA contribution of 35%, presumably because of much 

429 higher BC emissions in EA. The imported BC to EA comes mostly (65%) from SA, 

430 followed by those from EU (e.g., 30%). For the import of SO4+POM to EA, SA and EU 

43 1 emissions make quite comparable contributions, which is more than a factor of 3 larger 

432 than NA emissions. In SA, the imported SO4+POM is significant, with RAIR of 24%, of 

433 which 66% comes from EU, followed by 28% from EA. For BC aerosols in SA, the 

434 import from EA and EU emissions contributes more or less equally (e.g., 48% and 44%). 

435 

436 Our multi-model estimates of the relative contributions of local emissions and ICT import 

437 are generally consistent with results in literature. In NA, we estimate that local emissions 

438 account for about 81% of anthropogenic AOD, which is higher than 66% estimated by a 

439 modeling study [Leibensperger et al., 2012] . This difference would be narrowed if 

440 emissions in regions other than the four defined regions of this study are taken into 

441 account. Yu et al. [2012a] estimate based on satellite measurements that the imported 

442 pollution via trans-Pacific transport results in a clear-sky TOA direct forcing of -0.49 


19 



443 Wm" 2 over NA, which accounts for about 1 1.5% of total anthropogenic aerosol direct 

444 forcing. For comparison, our estimated contribution by ICT import to the direct forcing 

445 over NA is 14.7% in this study. The BC fractional contributions from different source 

446 regions as estimated in this study also agree well with a model estimate [Reddy and 

447 Boucher , 2009, referred to as RB09]. RB09 estimated that over North America, emissions 

448 from East and South Asia contribute to 18% and 6% of the total BC burden, respectively, 

449 which are somewhat smaller than our corresponding estimate of 23% and 8%. The 

450 contribution of North America emissions to BC burden in Europe was estimated at 5% in 

45 1 RB09, which is also somewhat smaller than our estimated 8%. RB09 estimated that the 

452 local BC emissions over South and East Asia accounts for more than 80% of the BC 

453 burden, which agrees well with our estimate of 82% and 84%. These comparisons show 

454 how the estimated significance of ICT import in this study agrees broadly with results in 

455 the literature, as these studies have defined source regions somewhat differently. 

456 

457 4. Conclusions and discussion 

458 We have assessed impacts of a 20% reduction of anthropogenic emissions in North 

459 America, Europe, East Asia, and South Asia on the aerosol optical depth and direct 

460 radiative forcing by using results from 10 global chemical transport models. On the basis 

461 of the multi-model average, a 20% reduction of anthropogenic emissions in the four 

462 regions combined lowers the global mean AOD and DRF by about 9%, 4%, and 10% for 

463 sulfate, organic matter, and black carbon aerosol, respectively. Despite the considerable 

464 model-to-model differences in the magnitude of forcing efficiency (direct radiative 

465 forcing per unit emissions), the models appear to consistently give a lower sulfate forcing 
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466 efficiency for SO 2 emissions from EA than that from EU, and the largest BC forcing 

467 efficiency for the EU emissions. These results may suggest that reducing BC emissions in 

468 EU would be relatively the most efficient way to mitigate the warming effect of BC 

469 aerosol, not considering cost and feasibility. However, a larger potential for BC 

470 reductions is available in East Asia. Reducing SO 2 emissions decreases the cooling effect 

47 1 of aerosols that counteracts the greenhouse wanning; but the impact of reducing EA SO 2 

472 is less pronounced (ca. 25%) than reducing EU, NA, or SA SO 2 by the same amount. 

473 However, such implications need to be further explored by taking into account a variety 

474 of aerosol impacts on climate through modifying cloud microphysics, changing 

475 atmospheric circulations, and altering the snow albedo. The air quality implication of 

476 any emission control should also be considered. 

477 

478 Our multi-model simulations of source-receptor relationships show that the impacts of 

479 the regional emission reductions are not confined to the region itself because of aerosol 

480 intercontinental transport. On an annual basis, intercontinental transport accounts for 10- 

481 30% of the overall AOD or DRF in a receptor region, compared to the influence of both 

482 regional emissions and intercontinental transport, depending on regions and species. 

483 While South Asia is most influenced by import of sulfate aerosol from Europe, North 

484 America is most influenced by import of black carbon aerosol from East Asia. Given that 

485 BC deposition may accelerate the melting of snow in the Sierra Nevada and cause water 

486 supply shortage in summer in the western U.S. [Hadley et al., 2010], the region may 

487 likely benefit from a future control of BC emissions in Asia. 

488 
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489 The multi-model assessment shows large differences between models in the impacts of 

490 emission reductions and the role of intercontinental transport, which highlights a need for 

49 1 improving models and developing observational databases for evaluating and 

492 constraining models. Quantifying anthropogenic AOD distributions from satellite 

493 measurements remains challenging, in particular over land. Current estimates of 

494 anthropogenic AOD based on total AOD and fine-model fraction measurements from 

495 MODIS are only feasible over ocean and subject to large uncertainties [Yu et al., 2009; 

496 Bellouin et al., 2005, 2008] . Such an estimate would be better constrained by a 

497 synergistic use of aerosol microphysical measurements as provided by other satellite 

498 sensors [Yu et al., 2012b and references therein]. Large-scale in situ measurements of 

499 anthropogenic AOD will be essential for evaluating satellite observations and model 

500 simulations. From the perspective of model improvements, efforts should focus on not 

501 only emission inventories [Textor et al., 2007] but also a variety of atmospheric 

502 processes that determine the atmospheric evolution of aerosols, such as parameterization 

503 of dry deposition and wet scavenging [Prospero et al, 2011]. Given that the aerosol 

504 intercontinental transport and its influences involve a wide span of scales, it is necessary 

505 to develop modeling systems that link the local, regional, intercontinental, and global 

506 scales. 
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705 Table Captions 

706 

707 Table 1. List of 10 models that participated in HTAP aerosol S/R experiments and are 

708 used in this analysis. 

709 

710 Table 2. Comparisons of AOD at 550 mn (x) from 9 HTAP models in this study with that 

711 inferred from Schulz et al. (2006). For SO4 and POM, AOD is inferred by using 

712 anthropogenic AOD and anthropogenic fraction of present day AOD. For BC, we infer 

713 AOD from the anthropogenic absorptive AOD reported in Schulz et al. [2006] by 

714 assuming that all BC is anthropogenic and BC has a single scattering albedo of 0.2 at 

715 550 mn. 

716 

717 Table 3. Annual mean aerosol direct radiative forcing normalized by AOD at 550 mn 

718 (NDRF, Wm'V 1 ) at top of atmosphere (TOA) and surface for SO 4 , POM, and BC) in 

719 the four source regions (NA, EU, EA, SA) and globe, which is derived from 

720 GOCART simulated monthly AOD and DRF for 2001. 

721 

722 Table 4. Change of global annual mean AOD (x x 1000, mean ± std. dev) in response to 

723 the 20% reduction of anthropogenic emissions in 4 source regions (SR6 - SRI) of 9 

724 HTAP models. 

725 

726 Table 5. Change of global annual mean aerosol DRF (unit: mWm' 2 , mean ± standard 

727 deviation, 1 mWm 2 = 0.001 Wm" 2 ) in response to the 20% reduction of anthropogenic 

728 emissions in 4 source regions as derived from analysis of SRI and SR6 runs of 9 

729 HTAP models. 

730 

73 1 Table 6. Global annual average TOA all-sky forcing efficiency relative to emissions 

732 from the source regions (unit: mWm' 2 per Tg) as derived from 8 HTAP models. The 

733 forcing efficiency for SO 4 and POM is calculated with respect to SO 2 and primary 

734 POM emissions, although some models with fully coupled chemistry include changes 

735 of SO4 resulting from the reductions of other emissions and a fraction of POM is 

736 secondary aerosol produced from a variety of volatile organic carbon. 
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737 

738 Table 7: Response of annual AOD and all-sky DRF (mean ± std. dev.) in the receptor regions 

739 (columns) to 20 % reductions of anthropogenic gas and aerosol emissions in the source regions 

740 (rows) estimated from 9 HTAP models. Percentage contributions of individual source regions 

741 are shown in parentheses for AOD (corresponding percentages for DRF are only slightly 

742 different and hence not shown). Contributions of a region on itself (termed “domestic”) are 

743 shown in bold. Here BC is listed separately from SO 4 and POM combined (in which SO 4 

744 makes a dominant contribution of 69-94%, depending on region), because of their distinct 

745 effects on solar radiation. For SO4+POM, surface DRF is slightly greater than TOA DRF 

746 because of weak absorption and not shown. A positive change of DRF for SO4+POM suggests 

747 a weakened cooling effect by the emission reduction. A BC emission reduction introduces a 

748 negative change of TOA DRF (a weakened warming effect) and a positive change of surface 

749 DRF (a weakened cooling effect). 

750 

751 
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752 

753 Figure Captions 

754 

755 Figure 1: Illustration of four HTAP defined regions for examining the source-receptor 

756 relationships for anthropogenic aerosols: North America (NA, 15°-55°N, 60°W- 

757 125°W), Europe (EU, 25°-65°N, 10°W-50°E), East Asia (EA, 15°-50°N, 95°-160°E), 

758 and South Asia (SA, 5°-35°N, 50°-95°E). Regional and annual anthropogenic 

759 emissions of SO 2 , primary POM, and BC from 8 HTAP models (excluding 

760 HADGEM2 and ECHAM5) are shown in bar charts, with error bar indicating the 

761 range of 8 models. 

762 

763 Figure 2: Annual average AOD for the external mixture of sulfate, POM, and BC 

764 simulated by baseline runs of 9 HTAP models. 

765 

766 Figure 3: Comparison of HTAP SRI 9-model median AOD for the external mixture of 

767 sulfate, POM, and BC (a) with MODIS-derived anthropogenic AOD over ocean as 

768 described in Yu et al. [2009] (b) and Bellouin et al. [2008] (c). 

769 

770 Figure 4: Comparisons of zonal variations of seasonal and 20°N-60°N average AOD for 

771 sulfate, POM, and BC combined as simulated by HTAP models (black line for median 

772 and shaded area for the range of 9 HTAP models) with the MODIS-derived over- 

773 ocean anthropogenic AOD (red line for Yu et al., 2009 and blue line for Bellouin et al., 

774 2008). 

775 

776 Figure 5: Seasonal variations of 9-model average TOA DRF (a), and surface DRF (b) in 

777 all-sky conditions (Win -2 ) for the external mixture of sulfate, POM, and BC as derived 

778 from the HTAP baseline simulations (SRI). 

779 

780 Figure 6: Global annual mean forcing efficiency (unit: mWm" 2 per Tg) for sulfate (a) and 

781 BC (b) with respect to regional (denoted as NA, EU, EA, and SA) emissions as 

782 simulated by 8 HTAP models. 
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783 

784 Figure 7: Annual average AOD (xlOO, top panel) and DRF (mWm' 2 ) at the top of 

785 atmosphere (TOA) (middle panel) and at the surface (bottom panel) in all-sky 

786 conditions resulting from 20% reductions of regional anthropogenic emissions over 

787 North America (NA), Europe (EU), East Asia (EA), and South Asia (SA), 

788 respectively. A positive value for DRF represents the reduced aerosol direct radiative 

789 forcing resulting from the reduction of emissions. Each of the four source regions is 

790 overlaid on corresponding SR6xx-SRl maps where xx represents NA, EU, EA, and 

791 SA. Individual aerosol components are assumed to be mixed externally. The results 

792 are based on 9 HTAP models. 

793 

794 Figure 8: Seasonal variations of changes in the all-sky DRF (mWm' 2 ) at surface (SR6- 

795 SRI) due to 20% reduction of regional anthropogenic emissions derived from 9 HTAP 

796 models. Each of the four source regions is overlaid on corresponding SR6xx-SRl 

797 maps where xx represents NA, EU, EA, and SA. 

798 

799 Figure 9: Ten-model (including ECHAM5) derived average relative annual 

800 intercontinental response (RAIR) for aerosol optical depth (AOD) in four receptor 

801 regions, by chemical component. Standard deviations, as indicated by error bars, 

802 reflect the model variability in simulating aerosol transport. 

803 

804 Figure 10: same as Figure 9 but for TOA all-sky DRF. 

805 
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Table 1. List of 10 models that participated in HTAP aerosol S/R experiments and are 
used in this analysis 

Model 

Model 

version 

Investigator(s) 

Resolution 
(Ion. x lat.) 

Major Reference(s) 

CAM-CHEM 

v3311ml3 

P. Hess 

2.5°xl.875° 

Pfister et al., 2008 

CAM-CHEM 

v3514 

ECHAM5 

HAMMOZ 

v21 

I. Bey, G. 
Forberth 

2.813°x2.813° 

Pozzoli et al., 2008a 
and 2008b 

GISS PUCCINI 

modelEaer 

D. Shindell 

5°x4° 

Koch et al., 2005; 
Shindell et al., 2006 

GMI 

v02a 

H. Bian 

2.5°x2° 

Bian et al., 2009 

GOCART 

v4p2 

M. Chin, T. 
Diehl 

2.5°x2° 

Chin et al., 2002, 
2007, 2009 

HADGEM2 

A-vOl 

N. Bellouin 

1.875°xl.25° 

Bellouin et al., 2011 

INCA 

v2 

M. Schulz 

3.75°x2.5° 

Schulz, 2007; 
Hauglustaine et al., 
2004; Textor et al., 
2006 

LLNL IMPACT 

T5a 

D. Bergmann, 
C. S. Atherton 

2.5°x2° 

Rotman et ah, 2004. 

SPRINT ARS 

v356 

T. Takemura 

1.125°xl.l25° 

Takemura et al., 2005 

Acronyms for model names: 

c Model - Chemistry version (NCAR, USA) 

ititute for Meteorology Hamburg Climate Model - version 5 with 

Kiel and MOZart chemistry sub-Module 

Space Studies, Physical Understanding of Composition-Climate 

acts model (NASA GISS, USA) 

A GSFC, USA) 

1 Radiation and Transport (NASA GSFC, USA) 

/ironment Model version 2 (Met Office, UK) 

;rosol (laboratoire des Sciences du Climat et de l'Environnement, 

Parallel Atmospheric Chemical Transport model (Lawrence 
laboratory, USA) 

sport Model for Aerosol Species (Kyushu University, Japan) 

CAM-CHEM: Community Atmospheri 
ECHAM5-HAMMOZ : Max-Planck Ins 
Hamburg Aerosol Me 
GISS PUCCINI: Goddard Institute for ! 

INteractions and Imp 
GMI: Global Modeling Initiative (NAS 
GOCART : Goddard Chemistry Aeroso 
HADGEM2: Hadley Centre Global Em 
INCA: Interaction of Chemistry and Ae 
France) 

LLNL IMPACT : Integrated Massively 
Livermore National I 
SPRINTARS: Spectral Radiation-Tran: 
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Table 2. Comparisons of AOD at 550 nm (x) from 9 HTAP 
models in this study with that inferred from Schulz et al. 
(2006). For SO 4 and POM, AOD is inferred by using 
anthropogenic AOD and anthropogenic fraction of present day 
AOD. For BC, we infer AOD from the anthropogenic 
absorptive AOD reported in Schulz et al. [2006] by assuming 
that all BC is anthropogenic and BC has a single scattering 
albedo of 0.2 at 550 nm. 



This study 

Schulz et al. [2006] 

tS04 

0.0335 ±0.0134 

0.0345 ±0.0164 

TpOM 

0.0108 ±0.0047 

0.0151 ±0.0094 

Tbc 

0.0022 ±0.0010 

0.0024 ±0.0010 
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Table 3. Annual mean aerosol direct radiative forcing normalized by AOD at 550 mn 
(NDRF, Wm'V) at top of atmosphere (TOA) and surface for SO 4 , POM, and BC) in the 
four source regions (NA, EU, EA, SA) and globe, which is derived from GOCART 
simulated monthly AOD and DRF for 2001. 

Sky 

Condition 

Region 

NDRFso 4 

NDRFpom 

NDRFbc 

TOA 

Surface 

TOA 

Surface 

TOA 

surface 

All Sky 

NA 

-24.9 

-24.5 

-28.4 

-39.1 

84.6 

-225.2 

EU 

- 21.1 

-20.4 

- 22.6 

-32.8 

93.0 

-190.8 

EA 

-21.4 

- 21.2 

-25.4 

-35.1 

83.2 

- 210.0 

SA 

-24.9 

-25.1 

-28.7 

-41.0 

89.6 

-235.2 

globe 

-24.2 

-24.1 

-30.0 

-41.5 

85.9 

-231.6 

Clear Sky 

NA 

-30.6 

-30.1 

-35.4 

-45.9 

62.1 

-246.1 

EU 

-25.6 

-24.6 

-27.9 

-37.9 

80.2 

-207.8 

EA 

-26.6 

-26.3 

-32.1 

-41.5 

59.9 

-231.3 

SA 

-27.7 

-27.8 

-32.7 

-44.8 

74.6 

-247.6 

globe 

-30.2 

-29.9 

-37.0 

-48.4 

62.8 

-252.7 
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Table 4. Change of global annual mean AOD (x x 1000, mean ± std. dev) in response to 
the 20% reduction of anthropogenic emissions in 4 source regions (SR6 - SRI) of 9 
HTAP models. 


Source Region 

^S04 

Tx>OM 

Tbc 

^SCM+POM+BC 

NA 

-0.67 ± 0.23 

-0.054 ± 0.036 

-0.026 ±0.013 

-0.75 ± 0.23 

EU 

-1.13 ±0.47 

-0.084 ± 0.058 

-0.050 ±0.021 

-1.28 ±0.45 

EA 

-0.97 ± 0.38 

-0.142 ±0.101 

-0.086 ± 0.047 

-1.19 ±0.37 

SA 

-0.32 ±0.11 

-0.102 ±0.059 

-0.038 ±0.019 

-0.46 ±0.10 

Total 

-3.10 ± 1.11 

-0.383 ± 0.246 

-0.199 ±0.093 

-3.70 ± 1.08 
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Table 5. Change of global annual mean aerosol DRF (unit: mWm' 2 , mean ± standard 
deviation, 1 mWm' 2 = 0.001 Wm~ 2 ) in response to the 20% reduction of anthropogenic 
emissions in 4 source regions as derived from analysis of SRI and SR6 runs of 9 HTAP 
models. 

Source Region 

DRFso4 

DRFpom 

DRFbc 

DRFs04+POM+BC 

All-sky TOA DRF (mWm' 2 ) 

NA 

15.4 ±5.1 

1.4 ±0.9 

-3.0 ± 1.6 

13.8 ±5.3 

EU 

24.4 ± 8.9 

1.6 ± 1.2 

-5.9 ±2.6 

20.2 ± 9.5 

EA 

19.4 ±7.9 

3.1 ±2.2 

-9.4 ±6.1 

13. 0± 9.9 

SA 

6.8 ± 1.9 

2.5 ± 1.7 

-3.6 ± 1.8 

5.7 ±2.6 

Total 

66.1 ±22.8 

8.6 ±5.7 

-21.9 ± 11.5 

52.8 ±25.3 

All-sky surface DRF (mWm 2 ) 

NA 

14.9 ±4.9 

1.9 ± 1.3 

5.5 ±2.4 

22. 4± 5.5 

EU 

23.0 ±8.0 

2.4 ± 1.7 

8.8 ±3.1 

34.1 ±7.3 

EA 

18.9 ±7.5 

4.4 ±3.2 

16.5 ± 8.7 

39.9 ± 11.7 

SA 

6.9 ±2.0 

3.6 ±2.3 

8.1 ±3.6 

18.6 ±4.6 

Total 

63.7 ±21.4 

12.3 ± 8.1 

39.0 ± 16.6 

115.0 ±26.0 
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Table 6. Global annual average TOA all-sky forcing efficiency relative to 
emissions from the source regions (unit: mWm ' 2 per Tg) as derived from 8 
HTAP models. The forcing efficiency for SO 4 and POM is calculated with 

respect to SO 2 and primary POM emissions 

, although some models with fully 

coupled chemistry include changes of SO 4 resulting from the reductions of 
other emissions and a fraction of POM is secondary aerosol produced from a 

variety of volatile organic carbon. 



Source Region 

S0 4 

POM 

BC 

NA 

-3.6±1.0 

-4.3± 1 .5 

27.2±14.2 

EU 

-3.7±1.1 

vq 

+l 

<N 

^t- 

36.7± 1 8.0 

EA 

-2.7±0.9 

-3.6±1.7 

27.7±19.0 

SA 

-3.8±1.0 

-4.0±1.7 

25.0±13.6 

Total 

-3.3±0.9 

-3.9±1.6 

29.0±16.8 


824 

825 


39 




826 


Table 7: 

Response of annual AOD and all-sky DRF (mean ± std. dev.) in the receptor regions 

(columns) to 20 % reductions of anthropogenic gas and aerosol emissions in the source regions 
(rows) estimated from 9 HTAP models. Percentage contributions of individual source regions are 
shown in parentheses for AOD (corresponding percentages for DRF are only slightly different 
and hence not shown). Contributions of a region on itself (termed “domestic”) are shown in bold. 
Here BC is listed separately from SO 4 and POM combined (in which SO 4 makes a dominant 
contribution of 69-94%, depending on region), because of their distinct effects on solar radiation. 
For SO 4 +POM, surface DRF is slightly greater than TOA DRF because of weak absorption and 
not shown. A positive change of DRF for SO 4 +POM suggests a weakened cooling effect by the 
emission reduction. A BC emission reduction introduces a negative change of TOA DRF (a 
weakened warming effect) and a positive change of surface DRF (a weakened cooling effect). 


Receptor Region 

Source 

Region 

NA 

EU 

EA 

SA 

SO 4 +POM AOD (xlOOO) 

NA 

-8.53i2.85 (82%) 

-1.22±0.66 (7%) 

-0.35i0.25 (2%) 

-0.24i0.22 (2%) 

EU 

-0.42±0.21 (4%) 

-13.81±4.43 (84%) 

-1.21i0.68 (6%) 

-2.18i0.79 (16%) 

EA 

-0.92±0.59 (9%) 

-0.65±0.48 (4%) 

-16.42i5.02 (84%) 

-1.05i0.66 ( 8 %) 

SA 

-0.54±0.91 (5%) 

-0.72±1.32 (4%) 

-1.68i0.95 (9%) 

-10.37i3.21 (75%) 

BC AOD (xlOOO) 

NA 

-0.28±0.09 (65%) 

-0.05±0.04 ( 8 %) 

-0.02i0.01 (1%) 

-0.01i0.01 (1%) 

EU 

-0.02±0.02 (4%) 

-0.51±0.15 (76%) 

-0.05i0.03 (4%) 

-0.08i0.04 ( 8 %) 

EA 

-0.10±0.09 (23%) 

-0.07i0.08 (10%) 

-1.04i0.28 (84%) 

-0.09i0.06 (9%) 

SA 

-0.04±0.03 ( 8 %) 

-0.04i0.03 ( 6 %) 

-0.14i0.08 (11%) 

-0.87i0.36 (82%) 

SO 4 +POM all-sky TOA DRF (mWnf 2 ) 

NA 

207.7±68.9 

28.8il4.9 

7.5i5.2 

6.3i5.7 

EU 

10.4±5.42 

304.4i93.3 

25.4il4.1 

54.2i20.3 

EA 

23.8±15.1 

14.4il0.9 

315.2i98.8 

25.0il5.5 

SA 

1 1.2± 15.5 

8.8i9.4 

32.9il5.4 

256.4i80.4 

BC all-sky TOA DRF (mWm‘ 2 ) 

NA 

-23.8±7.2 

-5.1i3.3 

-1.4il.3 

-1.5il.3 

EU 

-1.7±1.4 

-45.6il2.8 

-4.5i2.7 

-9.2i4.4 

EA 

-9.0±8.5 

-6.7±7.3 

-83.5i22.3 

-8.6i6.5 

SA 

-3.0±2.3 

-4.1i2.7 

-U.8i5.3 

-70.7i29.6 

BC all-sky surface DRF (mWm' 2 ) 

NA 

64.U20.5 

11.5±7.2 

3.3i3.0 

3.5i2.8 

EU 

4.3±3.5 

103.0i27.1 

10.0i5.9 

19.8i9.9 

EA 

22.7±21.6 

14.4il5.5 

212.8i57.0 

22.3il5.0 

SA 

7.9±6.2 

7.5±5.3 

29.3il3.5 

202.0i85.4 


40 




828 



r t / ^ 






NA 


EU 

A- 

w 








S0 2 (Tg/yr) 

50 i 

POM (Tg/yr) 

8 

2 

BC (Tg/yr) 

40 | 

30 

: 1 1 

1 . 

6 

4 

0 i 1 

i i 

1 

0 i 

ill 

L 

NA EU 

EA SA 

NA 

EU EA SA 

NA 

EU EA 

SA 


Figure 1: Illustration of four HTAP defined regions for examining the source-receptor 


relationships for anthropogenic aerosols: North America (NA, 15°-55°N, 60°W-125°W), 
Europe (EU, 25°-65°N, 10°W-50°E), East Asia (EA, 15°-50°N, 95°-160°E), and South 
Asia (SA, 5°-35°N, 50°-95°E). Regional and annual anthropogenic emissions of SO 2 , 
primary POM, and BC from 8 models (excluding HADGEM2 and ECHAM5) are shown 
in bar charts, with error bar indicating the range of 8 models. 
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Figure 2: Annual average AOD for the external mixture of sulfate, POM, and BC 


simulated by baseline runs of 9 HTAP models. 


42 





835 


836 

837 

838 

839 

840 

841 

842 

843 


(a) HTAP-median 

DJF 



JJA 



SON 





.00 .04 .08 .12 .16 .20, .30 .50 1.0 2.0 


(b) MODIS-YU09 
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Figure 3: Comparison of HTAP SRI 9-model median AOD for the external mixture of 
sulfate, POM, and BC (a) with MODIS-derived anthropogenic AOD over ocean as 
described in Yu et al. [ 2009 ] (b) and Bellouin et al. [ 2008 ] (c). 
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Figure 4: Comparisons of zonal variations of seasonal and 20°N-60°N average AOD for 
sulfate, POM, and BC combined as simulated by HTAP models (black line for median and 
shaded area for the range of 9 models) with the MODIS-derived over-ocean anthropogenic 
AOD (red line for Yu et al., 2009 and blue line for Bellouin et al., 2008). Note that HTAP 
model results cover both land and ocean, while MODIS-based anthropogenic AOD were 
estimated only over ocean. 
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(a) TOA DRF (b) Surface DRF 

DJF DJF 



-30 -20 -10 -8 -6 - 4 , -3 - 2-1 05 

Figure 5: Seasonal variations of 9-model average TOA DRF (a), and surface 
DRF (b) in all-sky conditions (Wm 2 ) for the external mixture of sulfate, POM, 
and BC as derived from the HTAP baseline simulations (SRI). 
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Figure 6: Global annual mean forcing efficiency (unit: mWnf 2 per Tg) for sulfate 
(a) and BC (b) with respect to regional (denoted as NA, EU, EA, and SA) 
emissions as simulated by 8 HTAP models. 


46 



852 


SR6NA-SR1 SR6EU-SR1 SR6EA-SR1 SR6SA-SR1 

AOD (xlOO) AOD (x 1 00) AOD (xlOO) AOD (xlOO) 



-20 -10 -8 -6 -4 -2 -1 -.5 -.2 -.1 0 AOD(xlOO) 


TOA RF (xlOOO) 



Surface RF (xlOOO) 



TOA RF (xlOOO) 



Tjt- A • 

• > 




Surface RF (xlOOO) 



TOA RF (xlOOO) 



Surface RF (xlOOO) 



TOA RF (xlOOO) 



§3 

Surface RF (xlOOO) 



0 5 10 20 50 100 200 400 600 1000 2000 5000 mlm" 

Figure 7: Annual average AOD (xlOO, top panel) and DRF (mWnT 2 ) at the top of 


atmosphere (TOA) (middle panel) and at the surface (bottom panel) in all-sky conditions 


resulting from 20% reductions of regional anthropogenic emissions over North America 


(NA), Europe (EU), East Asia (EA), and South Asia (SA), respectively. A positive value for 


DRF represents the reduced aerosol direct radiative forcing resulting from the reduction of 
emissions. Each of the four source regions is overlaid on corresponding SR6xx-SRl maps 
where xx represents NA, EU, EA, and SA. Individual aerosol components are assumed to be 
mixed externally. The results are based on 9 HTAP models. 
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Figure 8: Seasonal variations of changes in the all-sky DRF (mWm' 2 ) at surface (SR6- 
SR1) due to 20% reduction of regional anthropogenic emissions derived from 9 HTAP 
models. Each of the four source regions is overlaid on corresponding SR6xx-SRl maps 
where xx represents NA, EU, EA, and SA. 
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Figure 9: Ten-model (including ECHAM5) derived average relative annual 
intercontinental response (RAIR) for aerosol optical depth (AOD) in four receptor regions, 
by chemical component. Standard deviations, as indicated by error bars, reflect the model 
variability in simulating aerosol transport. 


856 

857 

858 


49 


859 


S04 


POM 


BC 


Mixture 


50% 


40% 

h 

PC 

O 30% 
u 

a 

£ 20 % 


10 % 


0 % 


T 

i 

! . 

T.lT 7 

..T TI 

J 1 

r nr 

1 II 1 

T 

,r l 

1 


NA 


EU 


EA 


SA 


Figure 10: same as Figure 9 but for TOA all-sky DRF. 

860 


50 


